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We extend a generic class of systems which have previously 

■ been shown to spontaneously develop scaling (power law) dis- 
^ | . tributions of their elementary degrees of freedom. 

Q-f While the previous systems were linear and exploded expo- 

00 . nentially for certain parameter ranges, the new systems fulfill 

y—{ , nonlinear time evolution equations similar to the ones encoun- 

1 — I ■ tered in Spontaneous Symmetry Breaking (SSB) dynamics and 

> : 

^NJ , evolve spontaneously towards "fixed trajectories" indexed by 

the average value of their degrees of freedom (which corresponds 

On i 

to the SSB order parameter). The "fixed trajectories" dynamics 
evolves on the edge between explosion and collapse/extinction. 

The systems present power laws with exponents which in a 
wide range (a < —2.) are universally determined by the ratio 

■ between the minimal and the average values of the degrees of 
freedom. The time fluctuations are governed by Levy distribu- 

> ■ 

tions of corresponding power. For exponents a > —2 there is 
no "thermodynamic limit" and the fluctuations are dominated 
by a few, largest degrees of freedom which leads to macroscopic 
fluctuations, chaos and bursts/intermitency. 
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The wide-spread appearance of scaling laws in various natural systems has attracted 
the attention of scientists for the last century [Q, ||. While for a measure zero 
set of systems (statistical mechanics systems in equilibrium at critical values of the 
parameters) scaling was thoroughly understood ||, ||, ||, for most of the scaling 
systems in nature, the origins of scaling remains a mystery J7|. It was felt for the last 
decade that in many cases the emergence || of scaling has to do with complexity and 
self organization, i.e. some systems have a tendency to organize spontaneously into 
critical states 0. Moreover, the identification of the relevant degrees of freedom of a 
complex system was shown to be essential in understanding its multiscale dynamics 



TO] in terms of its "irreducibly complex cores" |TT 



In a recent publication it was shown that the appearance of the scaling power 
laws is as generic in multiplicative stochastic systems as the Boltzmann law is in additive 
stochastic systems. By using the logarithm of the normalized degrees of freedom, we 
showed that even non-stationary systems can be investigated by mapping them into 
additive-stochastic systems similar to usual statistical-mechanics systems in thermal 
equilibrium. 

In the present paper we extend our mechanism to include a series of new features 
which makes it applicable to a wider range of phenomena. 

In particular we explain the natural emergence of "fixed trajectories" for which, (as 
opposed to the economics example in Jl2|) the average mean value does not diverge in 
time. 

These "fixed trajectories" are analog to the "fixed points" in equilibrium statistical 
mechanics systems and naturally evolve on the edge between exponential divergence 
and exponential decay. 

We show that on these "fixed trajectories", the systems display power laws and 
Levy fluctuations distributions which have wide universality basins with computable 
exponents. 
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For exponents a > —2 of the Levy distributions, one expects the inexistence of a 
"thermodynamic limit" at large system sizes and the appearance of bursts and intermi- 
tency dominated by the fluctuations of the values of a few individual degrees of freedom 
(low dimensional chaos). 

We call our generic mechanism Spontaneous Scaling Emergence (SSE) as it does 
not require criticality or self-organization to achieve the scaling state. The effective 
time evolution equations have common features with the equations of Spontaneous 
Symmetry Breaking. 

Before presenting the extensions let us discuss from the novel point of view certain 



aspects of the systems introduced in [13]. These systems are characterized by the 



following stochastic time evolution equation : 

w t (t + 1) = Xw t (t); i = 1,2,..., N (1) 
For defmiteness one can think of the example of economic systems [I2"[, where Wi(t) 



represents the wealth of investor i at time t. A is a random variable drawn from a 



generic distribution TI(A, t) . It was imposed in [|T^] that there is some lower bound to 
the value of Wi(t): 

Wi(t)>w w(t) Vi,t (2) 

where 

w(t) = l/Nj2Mt)- (3) 

i=l 

This dynamics was shown to lead to a power law distribution for the various w values : 

$(w,t) = C(t)w a (4) 

where C(t) is a normalization factor and a = —1 — ^ . For values of Wq small 
compared to 1/N (i.e. — log^o > logiV) there are corrections to (f|) originating in the 
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contributions of large w values. The large w corrections are relevant only for power 



laws with exponents a more positive than —2 and were discussed in []T3| |. For the sake 
of brevity we will not discuss them here again in the new context. One also expects a 
deviation from (Q) in the immediate vicinity of the lower cut-off wq w(t). 

With these restrictions, the striking property of (Q) is its independence on n(A, t) 
except for the coefficient C (t) which might depend on time in a quite irregular fashion. 

Since we are looking here for systems in which C (t) flows naturally and robustly to 
finite C values, it is instructive to recall in some detail how the robustness of arrised 
in The crucial step was to consider in place of the time-dependent distribution 
(|j) the distribution of the normalized variables 

uj{t) = w t /w(t) (5) 

The distribution P(u,t) of the variables u(t) is "pinned" at least against overall 
rescalings to — > K{t)uo by the condition 

J ' ujP(uj,t)duj = 1 (6) 

which follows trivially from the normalization ([5]). Equation ([!]) becomes in terms of 
the new variables: 

u i (t + l) = fiuj i (t) (7) 

where the distribution H'([i) of the random variable fi is related to the distribution 
n(A) of A through the relation 

^// \ , u)(t + 1) . . 

This overall rescaling in the argument of II turned out to be enough to eliminate the 
dependence of C on t in (|j). This means that the use of uj's in place of w's rescales in 
a very particular way the argument of II. More precisely U'(fj,,t) is fine-tuned by the 
transformation (|5|) such that P(u,t + 1) = P(u,t), i.e. (eq. 9 in Ref [12]) 



|n / (/i,t)/i 1+1 /( 1 -^rf / i = 1 (9) 

This steady state situation is obviously reached by having the distribution of the 
fi's which appear in (0) finely tuned to equilibrate the values of \i larger than 1 (which 
lead to the increase of the a>'s) with the smaller /i values (which decrease the u'a). 

One sees now that one could in principle allow the shape of the distribution n(A, t) 
to depend on w as long as an overall rescaling of its parameter A — > o(t+i) ^ i nsures that 
the distribution stability condition (||) holds for each u i.e. / n'(/i, t, u;)/i 1+1// ( 1_w °)(i// = 
1 Vtu, t. In terms of the stock market this means that the power law is stable even if 
the spread (width) of the gain-loss distribution depends on the wealth of the investor 
(or varies in time). All that is required is that the average gain expectancies of all the 
investors are equal (except for the poorest, which can not lose more than allowed by 
the poverty bound wq). 

This is quite a large degree of universality even though unlike critical phenomena 
of statistical mechanics the scaling exponents depend here on a continuous parameter. 

Now it is clear what is needed in order to get rid of the time dependence of C(t) in 
(|j): one needs to find a mechanism which brings n(A) to be finely balanced between 
A > 1 and A < 1 values, such as to fulfill equation (|S[). Such a mechanism will 
prevent the exponential explosion or collapse of (H). In the stock market example 
this mechanism is expressed by the fact that as long as the total real wealth of the 
community is unchanged, overall gain factors across the entire population are irrelevant 
(they express just inflation in the instruments of payment). It will turn out that this 
is a degenerate example of a wider generic phenomenon. 

In order to see how our new models exploit this idea, consider a system governed 
by the stochastic evolution rule : 

Wi(t + 1) -Ui(t) = K(wi(t) +u; ); * = 1,2, ...,N (10) 
£*;<(*) >0 Vi,t (11) 
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where k is a random variable drawn from a distribution T(k, u)(i)). We use different 
symbols for the analogues of A and II of ([!]) in order to avoid possible confusion. More- 
over note that in the new models there is no normalization condition on the elementary 
variables tOiif) and therefore their average u(t) is a variable. In fact, T(k, u)(t)) de- 
pends on u(t) and is assumed that for u(t) = the distribution T(k, 0) favors the 
appearance of strongly positive values of k while for u(t) at a larger scale W, the dis- 
tribution T(k,W) favors negative values of k. According to (|10D, this means that for 
small values of u)(t) the uj's will have the tendency to increase while for u(t) ~ W they 
will have a tendency to decrease. This means that there exists a value uo eq (determined 
by a condition of the type (^) towards which u(t) will always flow, by the virtue of 
(|l(J).(We defer to further publications the discussion of non-monotonic T's which lead 
to multiple solutions uj eq and to possible branchings/bifurcations). 

Once uj eq is reached, eq. (|IDD becomes of the type (0) (for the variable U)i(t) + u ). 
Independently of the details of T(/c, u>(t)) one then has : 

P{u) = C(u + Wo )-i-i/(i-^/o«) (12) 
with constant C. This is a stationary power law, similar in form to (f|), but with 



no time dependence ! According to the analysis in |12| , the time fluctuations of the 
system are then characterized by a random walk with steps of variable stochastic length 
distributed according to flT2|). This means , according to the generalized central limit 



theorem, that the time fluctuations of various quantities are given by a Levy distribution 



17| of corresponding power. 

In the case in which the Levy distribution is wider (and the the power law exponent 
a is > —2. due to the effect of large Wi values Hl3|), large macroscopic fluctuations 
are expected due to the variations of the large io/s. In particular, they might bring 
a system macroscopically away from the fixed point u eq . This leads to a burst of 
non-equilibrium activity until enough correcting (relaxation) steps take place. In the 
case that it exists more then one uj eq value, these macroscopic fluctuations may induce 
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jumps between the various "equilibrium" trajectories. The fact that for a > —2 Levy 
distributions only the largest elements are relevant to the dynamics predicts the absence 
of a thermodynamic limit as the number of elementary degrees of freedom is taken to 



infinity. Such a phenomenon might explain some of the simulation results in [14], [15]. 
We propose to study these effects numerically and compare with the theoretical and 
experimental universal properties of chaotic systems [16| . 

The assumption in fllCf ) that the distribution T(k, Lj(t)) depends on u(t) might seem 
arbitrary, but in fact it arises naturally in quite familiar contexts. Consider a generic 
complex system in which some cooperative "short-range" underlying dynamics insures 
that certain subsystems act coherently and can be represented by a single variable 
uJi(t). One can think of such subsystems as families, clans, populations within an 
ecology, etc. Assume that isolated from each other these systems act according to ([!]). 
Suppose however that there is an "infinite-range" competitive interaction between the 
u's (clans competing for land, populations for food, etc.). The resulting dynamics may 
then be of the form 

u}i{t + 1) — uJi(t) = KUi(t) + uj u(t) — -uj(t)uJi(t) + etc. , (13) 

where etc. stands for possible additional inhomogeneous terms. One sees that for 



small values of ui the stochastic equation (|T3|) drives the a>'s to higher values while for 



large values of u) the k/s are driven to smaller values and one is in the case fllCf ). 

Note that equation ( Jl3| ) is very similar to the familiar equation governing the emer- 
gence of Spontaneous Symmetry Breaking (SSB): 

M kM - ^M 2 + etc. (14) 

from which it differs only by the fact that it mixes the "short-range" variable uJi(t) with 
the "infinite-range" variable u(t). 

This suggests the extension of our models for "vectorial" u variables which take 
values in spaces with various symmetries. 



In terms of the SSB vocabulary, the power of the distribution flT2|) is governed by 
the ratio between the scale of the global variable u)(t) and the scale of the local variable 
uio- As seen above, both of them may be dynamically fixed by the interplay of various 
non- homogeneous terms in (|13|). Note that the emergence of power laws and hierarchies 
follows then naturally even if the ratio of the 2 scales is of order 1 (e.g. 0.2 — 0.3 in the 
economics application |L2|). 

In further studies one can consider the effects of higher powers in ( ]T3| ) and relax the 
"infinite-range" assumption. In fact, one can study multiscale hierarchies of ranges, and 
consider the domains of almost constant cj's as the elementary degrees of freedom Q of 
a coarser system. For instance, each family within a clan acts as a unit but competes 
with the other families for the clans resources. Yet at a larger scale, the clan acts as 
a self-reinforcing unit in competition with other clans over the country resources, and 



so on. Examples in diluted finite connectivity spin glasses [I I] and Weiss domains in 
magnetic systems may be of relevance. 

This multiscale cooperation-competition interplay might be particularly relevant for 
the study of the speciation process in various ecological and territorial conditions. It 
also raises some suggestions on the ways to implement the reductionist program as 
one passes from one scale to another. The scaling laws and their relevant degrees of 
freedom might give some objective criterion for identifying the passage from one level 



of description to another. In the spirit of [10| and | L| one can "turn the tables" and use 



the power laws in order to identify from (|T2"|) the size and identity of the relevant effective 
degrees of freedom of the system: the "irreducibly complex cores" in the terminology 
of jllj]. In turn, the dynamics of uj would indicate the "multiscale-complex" dynamics 
of the hierarchy of macros in the terminology of [[H] . 



Let us reiterate in conclusion the main points discussed in the present paper. 

1. We have shown that systems with local self-catalizing multiplicative interactions 
and competitive global interactions (|10"D (|13D , may generically flow into "fixed 
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trajectories" at the edge between explosion and collapse. 

2. The various "fixed trajectories" are indexed by the values of the average degrees 
of freedom which emerge from a nonlinear time evolution equation similar in form 
to the one governing SSB. 

3. In those equations, both local and global variables appear and their scales (uq 
and Co) may be fixed by the generic nonlinear dynamics. 



4. The degrees of freedom are distributed according power laws (|I2) whose exponents 



are determined universally (for a wide range) only by the ratio of the 2 scales ujq/uj. 

5. The fluctuations around the "fixed trajectories" are governed by Levy distribu- 
tions of corresponding (universal) exponents. 

6. In the cases in which the Levy distribution has exponent a > —2 the system does 
not have a thermodynamic limit and presents bursts and intermitency produced 
by the few degrees of freedom with the largest values (chaos [0]). 

We are now facing the enormous task of applying the present framework to particu- 



lar complex systems beyond the economic example presented in ||12|| . The most urgent 
objective is to relate for each system the complex dynamics to the appropriate tu's, 
identify uo and relate it by ([12]) to the actual power law measured experimentally. 

The values of Uq may be determined by different mechanisms in various systems. 
Examples for such mechanisms are the extinction of species with too few members, 
subsidies to the poor, minimal size of population needed to found a city, discretization, 
additive noise, pumping of energy, etc. 
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Note Added in Proof 



A detailed discussion of some generalizations of the model introduced in [I^| can 
be found in "Convergent Multiplicative Processes Repelled from Zero" by D. Sornette 
and R. Cont. We thank D. Sornette for very instructive correspondence. We thank 
D. Stauffer for advice and encouragement since the beginning of our research on the 
present subject. 
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